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ABSTRACT 

We present Spitzer images of the Taurus Complex (TC). We take advantage 
of the sensitivity and the spatial resolution of the observations to characterize the 
diffuse infrared emission across the cloud. This work highlights evidence of dust 
evolution within the translucent sections of the archetype reference for studies of 
quiescent molecular clouds. We combine the Spitzer 160 /xm and IRAS 100 /xm 
observations to produce a dust temperature map and a far-IR dust opacity map 
at 5' resolution. The average dust temperature is about 14. 5K with a dispersion 
of ±1K across the cloud. The far-IR dust opacity is tightly correlated with the 
extinction derived from 2MASS stellar colors and is a factor 2 larger than the 
average value for the diffuse ISM. This opacity increase and the attenuation of 
the radiation field both contribute to account for the lower emission temperature 
of the large grains. The structure of the TC significantly changes in the mid-IR 
images that trace emission from PAHs and very small grains (VSGs). We focus 
our analysis of the mid-IR emission to a range of ecliptic latitudes away from the 
zodiacal bands and where the zodiacal hght residuals are small. Within this cloud 
area, there are no 8 and 24 fira counterparts to the brightest 160 fim emission 
features. Conversely, the 8 and 24 /im images reveal filamentary structure that is 
strikingly inconspicuous in the 160 /im and extinction maps. The infrared colors 
vary over sub-parsec distances across this filamentary structure. We compare the 
observed colors with model calculations quantifying the impact of the radiation 
field intensity and the abundance of stochastically heated particles on the dust 
SED. To match the range of observed colors, we have to invoke variations by a 
factor of a few of both the interstellar radiation field and the abundance of PAHs 
and VSGs. We conclude that within this filamentary structure a significant 
fraction of the dust mass cycles in and out the small size end of the dust size 
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distribution. 



Subject headings: 
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Introduction 



Understanding interstellar dust evolution is a major challenge underlying the 
interpretation of many infrared (IR) observations. The composition of interstellar dust is 
the outcome of processes in interstellar space that break and r e-build grains over tirae sca les 



much shorter than that of the dust renewal by stellar ejecta (jPraine 



Jones 



2OO5I : 



1990 



Tielens 



1998 



Zhukovska fc Gail 



20081 ). There is a wide consensus on such a conclusion, but 
the processes that drive dust evolution in space are still poorly understood. Understanding 
the processes that control the small grain end of the dust size distribution is of particular 
importance because the smallest dust particles have a marked impact on the physics of 
the ISM. Polycyclic Aromatic Hydrocarbons (PAHs), and carbon nanometer Very Small 
Grains (VSGs) absorb a significant fraction of the UV starlight and are considered to be the 
dominant gas heating source in interstellar space penetrated by UV light. Given their large 
total surface area, it is also proposed that they a ct on gas dep letion, chemistry, ionization 



degree and coupling with the magnetic field (e.g. iTielens 



20081). 



It is in molecular clouds where dust may grow by accretion of gas phase elements 
and grain-grain coagulation. Observations of the translucent sections {Ay less than a 
few magnitudes) of molecular clouds are of interest to trace dust evolution that occurs 
before matters enter dense molecular cores and form stars. Many UV and IR studies 
provide evidence of diverse processes that impact the dust size distribution. The systematic 
chang e in the visible-to-U V extinction curve correlated with the Ry = Ay/ EjB — V) 



ratio (ICardelli et al. 



19891 ) has been interpreted as evidence for grain growth (IKim et al. 



19941 ). Variations of the dust IR emission colors in translucent molecular clouds have 
been interpreted as evidence for changes in the small grains abundances and large grai ns 



emissivities (e.g. 



2005 



Kiss et al. 



Bernard et al. 



20061 ). 



1993 



1999 



Stepnik et al. 



del Burgo et al. 



2003 



Schnee fc Goodman 



2OO3I ) have shown that the deficit of small grains 



- 5 - 



and enhancement of large grain emissivity within a quiescent dense filament can be 
interpreted as grain-grain coagulation. However, some UV observations do not fit with 
a simple evolutionary picture where the smallest grains progressively stick on the largest 



( jBoulanger et al. 



1994 



Mathis 



1994 



Whittet et al. 



20041 ). 



Corr elations with density and 



example, 



Miville-Deschenes et al 



velocity structure of clouds have also been measured. For 



(I2OO2I ) observed an enhancement of PAHs abundance in a 



f ilame nt of Ursa Major characterized by a large transverse velocity gradient. 



Bernard et al 



(119991 ) measured a drop in the PAHs abundance between H I gas and denser molecular 
gas traced by CO emission. Mid-IR spectroscopic variations of photon-dominated regions 
suggest that photo-chemical processing of VSGs gives birth to free PAHs in PDRs 



(jRapacioli et al. 



2005; 



Berne et al 



2007 : 



Compiegne et al 



20081 ) 



Most of the earlier work focuses on sm a. 
low resolution data (e.g. iLagache et al.lll998 



1 sections of molecular clouds or is based on 



Cambresy et al. 



2OO5I ). In this paper, we 



take advantage of the coverage, sensitivity and resolution of recent Spitzer observations to 
present a global stud y of the nearby Taur us Complex (TC), a giant molecular complex about 



140 parsecs from us (IKenyon et al. 



19941 ) characterized by a lack of young massive stars and 



clusters. The TC is perfectly suited for high spatial resolution stud i es, mainly focused on 



1999 



Sunada fc Kitamura 



1999 



low mass star formation (e.g. iHartmann 
and interstellar chemistry (e.g 



2002 



'ratap et al. 



Harju et al 



Palla 



1997 



2000 



fc Stahlei 



2OO2I: 



Nakam 



Wolkovitch et al 



Whittet et al.l 



1997 



2007 



ura &: Li 



2008h 



Maezawa et al. 



Goldsmith et al 



20081 ). The TC also is the ideal region to analyze the structure of the interstellar medium 
(ISM) since it corresponds to the intermediary phase between diffuse, turbulent, non 
structured medium and the massive star forming regions dominated by the radiation field 
of young OB stars. Howe ver, only a few glob al stu dies have been perfo rmed through 



infrared (IR) observations. 



Langer et al 



fll989h and 



Abergel et al 



(I1994I ) have shown that 
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the far-IR (FIR) emission as seen by IRAS was correlated to ^^CO emission and extinction, 
and that the cold dust component (low 60 am to 100 fim emission ratio) was particularly 



well traced by the gas. 



Stepnik et al 



(120031 ) have observed significant variations of the IR 



colors within a TC filament which may be related to an evolution of the dust properties as 



one enters or exits the most exti nguis 



the TC by 



Whittet et al. 



led regions o 



(12004h and 



Shenoy et al. 



the c 



oud. Analysis of extinction within 



(120081 ) have also highlighted variations 



of dust properties. 

We present the Spitzer observations of the cloud and the way we process them to 
build complete maps of the extended emission in section [21 focusing on the zodiacal light 
subtraction issue. Section [3] presents our results on the average dust emission spectrum, the 
dust temperature variations and the correlation between FIR emission, dust temperature 
and visual extinction. We analyze a small sub-area of the TC in section H] to further 
constrain the local variations of the dust emission colors and relate them to variations of 
the ISRF intensity and dust size distribution in section [51 



Observations 



2.1. Spitzer observations 



The bulk of our analysis relies on recent data obtained by 



Spitzer using the I RAC and 



MIPS instruments and are part of the Taurus Legacy project ([Padgett et al.l 120091 ). Briefiy, 
the data come from three Spitzer programs (PIDs 3584, 30816 & 462), use two epochs to 
find and weed out asteroids, and cover a total area of approximately 44 square degrees. The 
IRAC data in their four bands (3.6, 4. 5, 5.8 fc 8 um) ut ilize the 12 seconds integration High 



Dynamic Range (HDR) frame mode (iFazio et al 



2004] ) with 2 dither positions per epoch. 



In the present study, we only use the IRAC 8 fim channel. The MIPS data in their three 
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Fig. 1. — MIR emission of the TC as seen by (a) IRAS 12 and (b) IRAC 8 /im at the 
resolution of IRAS 100 (4.3') and after point source subtraction. The units are in MJy/sr. 

bands (24, 70 & 160 /xm) utilize the fast scan mode with a 3 seconds integrat ion time per 



fram e and 5 visits per pixel. Because of the different size of the three arrays (iRieke et al 
20041 ) only the 24 /im maps have full coverage twice (two epochs). The 70 and 160 yum 
maps do cover the entire survey region, but only to a depth of 15 seconds integration per 
pixel. We use the three MIPS channels in the present study. 

As part of the Legacy commitment to improve the data, most of the obvious artifacts 
(latencies, stray light, cosmic rays, etc) for both IRAC and MIPS are removed in the 
individual frames (the Basic Calibration Data or BCD), before the creation of the final 
mosaics. Two additional steps are carried out also at the BCD level prior to mosaicking: 
zodiacal light removal and overlap co rrection. The f irst st ep uses the Spitzer Science Center 



(SSC) zodiacal light model based on 



Kelsall et al. 



(Il998l ). while the second step matches 



(globally) the median sky level across frames by applying a small offset (either positive or 
negative). Although these two steps rely on standard techniques and our best knowledge 
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Fig. 2.— MIR emission of the TC as seen by (a) IRAS 25 and (b) MIPS 24 /im at the 
resolution of IRAS 100 (4.3') and after point source subtraction. A zodiacal band appears 
clearly on both images at a different position due to the change in parallax. On the MIPS 24 
Hm images, offsets due to the combination of different epochs observations are also visible. 
The units are in MJy/sr. 



of the sky at mid/far infrared wavelengths, we do check processing does not modify the 
overall levels of the diffuse emission. Therefore, we compare ou r mosaics to the IRIS data 



(reprocessed IRAS data by 
(see Fig m El [3 and H. 



Miville-Deschenes fc Lagache 



20051 ) at 12, 25, 60 and lOO/im 



2.2. Zodiacal light residuals 



The TC is located at a low ecliptic latitude (/? < 10°). We use the model of iKelsall et al, 

(119981 ) to estimate the zodiacal light surface brightness. They significantly change in the 
multi-epoch Spitzer data. We only report the values of the first epoch: ~10 MJy/sr at 8 
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(a) 



(b) 



Fig. 3.— FIR emission of the TC as seen by (a) IRAS 60 and (b) MIPS 70 /im at the 
resolution of IRAS 100 (4.3'). The units are in MJy/sr. 
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Fig. 4.— FIR emission of the TC as seen by (a) IRAS 100 and (b) MIPS 160 /im at the 
resolution of IRAS 100 (4.3'). The units are in MJy/sr. 
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/iin, ~50 MJy/sr at 24 /iin, ~15 MJy/sr at 70 fim and ~5 MJy/sr at 160 fim. The total 
surface brightness observed by Spitzer are about 15, 55, 25 and 135 MJy/sr at 8, 24, 70 
and 160 fim respectively. The contribution of the zodiacal light to the infrared emission 
is significant at all wavelengths shortward of 160 /im. The subtraction of this foreground 
component is thus important but difficult, as discussed below. 

Subtraction of the zodiacal light component for the IRAC data is done as follows. 
Each 5 by 5 arcmin IRAC BCD image has an estimate of the zodiacal light a ssociated to 
i ts hea der (ZODY_EST keyword). This estimate is also based on the model of 



Kelsall et al. 



( 119981 ). However, the Spitzer pipeline subtracts from each BCD the contribution of the 
sky (including zodiacal light) at the sky bias observation position ( SKYDKRZB k eywor d). 



Kelsall et al. 



(119981) 



In order to fully subtract the zodiacal light in agreement with the 
model, we have to subtract the difference (ZODY_EST - SKYDKRZB) for each single BCD 
when building the mosaic. A correction factor for extended emission is also required, as 
given in the IRAC Data Handboolcl- Comparison with IRIS 12 /im images shows that 
some structures have not moved in more than 20 y ears (see Fig. HI) . We thus believe these 



Kelsall et al. 



mm . 



structures are not associated with zodiaca l band s (ILow et al. 
into account by the model of 



19841 ). which are not taken 



For MIPS 24 and 70 fim data the steps described above for IRAC are not necessary 
since the image bias is estimated a different way. The mosaic is built with a correction that 
takes into account the predicted value of the zodiacal light for each BCD. As for IRAC 8 
/xm, we compare MIPS 24 and 70 fim to IRIS 25 and 60 fim images (see Fig. [2] and [3]) . At 24 
and 25 fim, there are some structures that have not moved between the time they have been 
observed by IRAS and by Spitzer. However, it appears on the MIPS 24 /xm mosaic that 
there is an emission component of constant ecliptic latitude 2° < (3 < 3°. This component 



^See http:/ /ssc. spitzer. caltech.edu/irac/dh 
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is also visible on IRIS 25 fim images, but at a slightly different position, for 1° < (3 < 2° 
(see Fig. [2]). The positional difference of the features between epochs suggests that it is 
a zodiacal dust band whose position on the sky has changed between IRAS and Spitzer 
observations due to parallax. On the MIPS 24 /im mosaic, there are also some offsets that 
appear between the different epochs when the TC observations have been scheduled. We 
assume these offsets come from the uncertainty of the zodiacal light model. At 60 and 70 
fim, the emission structure is almost the same and we are confident that the contribution 
of the zodiacal light residuals, including bands, is neghgible. 

We estimate the residuals associated with the zodiacal band by median filtering the 
MIPS 24 fim mosaic as a function of the ecliptic latitude. The standard deviation of the 

R colors within the 



Kelsall et al. 



fll998h . 



resulting profile is about 0.3 MJy/sr. We compute the zodiacal light 
TC at the first two epochs of the observations thanks to the model of 
We use them to scale the MIPS 24 /im measurement to the other wavelengths. We obtain 
standard deviations of less than 0.03 MJy/sr for MIPS 160 fim, about 0.1 MJy/sr for MIPS 
70 fim and up to 0.5 MJy/sr for IRAC 8 fim. The spatial variations of the residual maps are 
about 21, 3.4, 0.47 and 0.15 MJy/sr at 160, 70, 24 and 8 fim respectively. The contribution 
of the zodiacal band is not significant at MIPS 160 fim or at MIPS 70 fim. However, it 
is an important component of the emission variations for MIPS 24 fim and IRAC 8 fim. 
Thus, while we will conduct the analysis of the extended emission across the entire TC at 
wavelengths strictly longer than 24 fim, we will adapt our analysis at shorter wavelengths 
in such a way we avoid the main zodiacal band. 



2.3. Point source subtraction 

We are interested in the extended emission across the TC. Therefore, we have to 
remove the point sources, which are mostly present at shorter wavelengths (IRAC 8 fim 
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and MIPS 24 /iin). We first produce a filtered map resulting from a median filter of the 
data by a box of 3'x3'. We subtract this median filtered image from the original mosaics. 
Finally, every pixel in the difference which is more than 3a away from the median value 
within a box of 3'x3'is flagged and replaced by interpolated value of neighboring pixels. 

The source removed Spitzer data at 8, 24, 70 and 160 /im are shown at the resolution 
of IRAS 100 /im (4.3') observations on Figures [T|, [2], [3] and HI For sake of comparison, IRIS 
images at 12, 25, 60 and 100 /xm are also shown. 



3. Data analysis 

We begin our analysis by comparing the infrared observations with each other on a 
global scale. Some of these observations can be related to one dust component, especially 
for Spitzer observations: IRAC 8 /im traces PAHs, MIPS 24 fim traces VSGs and MIPS 160 
/im traces BGs. The comparison between the shorter, intermediate and longer wavelengths 
observations provide us with informations regarding the variations of the PAHS, VSGs and 
BGs relative abundances, their properties and/or their excitation conditions. We combine 
Spitzer, IRIS as well as the Zodi-Subtracted Mission Average Maps of DIRBE to build an 
average sp ectral energy distribution (SED) of the TC and compare it to a diffuse medium 



reference ( iFlagey 



20071 ). We use this average SED to constrain average dust properties 
across the TC. We then analyze of spatial variations by building a map of the big grains 
temperature and FIR emissivity that is compared to FIR emission and visual extinction. 



3.1. Global comparison of infrared observations from 8 to 240 /xm 



We examine the available infrared observations from 8 to 240 /^m: IRAC 8 /im, MIPS 
24 /im, IRIS and DIRBE 60 /im, MIPS 70 /im, IRIS and DIRBE 100 /im, DIRBE 140 /im, 
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MIPS 160 /iin and DIRBE 240 /iin. The data (except DIRBE channels) are shown at the 
resolution of IRAS 100 fim observations (4.3') on Figures [T] to |H The emission morphology 
is significantly different between shorter and longer wavelengths with the 60, 70, 100 and 
160 fim data showing a different correlation and the 8 and 24 /im emission exhibiting a 
strong correlation. It thus appears that, away from the regions highly contaminated by the 
zodiacal bands, emission from PAHs and VSGs are well correlated with each other while 
they are not with that of BGs as further discussed in Section |H 



3.2. Global dust spectral energy distribution 



To build the SED, we plot the fiux at each wavelength as a function of the immediate 
longer wavelength (e.g. MIPS 70 fim versus IRIS 100 /im). Whenever possible, we measure 
the average color of the dust emission by fitting this pixel-to-pixel correlation by a straight 
line. Three separate IRIS plates cover th e entire TC and the correlations ar e done on a 



plate by plate basis to avoid offset issues (IMiville-Deschenes fc Lagache 



20051 ). In order to 



constrain efficiently the dust properties, we extend the SED as far as possible towards the 
longer wavelengths and add the DIRBE colors to Spitzer and IRIS measurements. 

Each plate exhibits a strong correlation between IRIS 60 fim and MIPS 70 fim, MIPS 
70 yum and IRIS 100 yum, IRIS 100 fim and MIPS 160 fim, which are all dominated by 
the BG component. The average dust emission colors are MIPS70/IRIS60 = 2.2 ± 0.3, 
IRIS100/MIPS70 = 2.9 ± 0.3, MIPS160/IRIS100 = 4.0 ± 0.15. The uncertainties on 
the TC colors take into account the variations from one IRIS plate to another. We also 
measure, at the resolution of DIRBE, the DIRBE60/DIRBE100, DIRBE100/DIRBE140 
and DIRBE140/DIRBE240 colors as well as the DIRBElOO/IRISlOO color to match DIRBE 
colors with the previous measurements (see Tab. [1]). We sh ow the a verag e FIR SED of the 



TC in Figure [5] along with that of the diffuse medium from 



FlagevI (120071 ) 
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Fig. 5. — Far infrared mean spectral energy distribution of the TC (stars) and the dust 
model that fits it [solid line) as well as the v'^By{T) function that best fits the FIR section 
of the SED. For the sake of comparison, we show the SED of the extended mean Galactic 
diffuse emission [diamonds) and the dust model that fits it [dashed). The TC average SED 
is normalized at the DIRBE 100 /im diffuse brightness. 

3.3. Dust temperature 



We use the average SED of the TC to constrain the BGs temperature, focusing on 
the wavelengths from 100 to 240 /im, where this component is expected to dominate the 
emission. We fit this reduced SED by a u^B„(T) function. We assume /? = 2, which 



is the value observ ed in the diffuse IS M (IBoulanger et al. 



large grains i n the 



Desert et al 



19961 ). It is also that of the 



(1l990l ) model and of graphite and silicate grains in the 
Draine fc Lil (120071 ) model. The best fit, shown on Figure [5] as a dotted line, corresponds 
to an average BGs temperature of 14.5 ± 0.2 K. This value is in agreement with previous 
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Fig. 6. — Dust temperature map as deduced from the IRIS 100 fim to MIPS 160 fim bright- 
nesses ratio and a v'^By{T) function. The contours are those from the MIPS 160 /im map 
smoothed to the IRAS 100 /im resolution of 4.3'. The slight level jumps in the map (e.g. at 
a declination of ~ 26.5°) are due to the offsets between IRIS plates. 



Stepnik et al. 



fl2003h 



ope temperature 



measurement of the BGs temperature within a filament of the TC by 
They have measured a large scale temperature of 16.8 K and an enve 
of 14 .8 K using IRAS 100 /xm and PRONAOS 580 /xm observations fiRistorcelh et al 
19981 ). They also noted a drop of the 60 /im emission relative to the 100 /xm that they 
explained by a depletion of the VSGs and by the introduction of a cold BGs component in 
addition to the warrn e r EG component. The cold component has a temperature of 12 K. 



del Burgo fc Laureijd (120051 ) used ISOPHOT and IRAS measurements from 60 to 200 /xm 
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Fig. 7. — Correlation plot of the dust temperature as deduced from IRIS 100 and MIPS 160 
/im observations as a function of MIPS 160 fim brightness. The plot corresponds to one IRIS 
plate only. 

to separate the dust emission into a cold and a warm component for the Taurus Molecular 
Cloud TMC-2 and its surrounding. Their temperatures are 12.5 and 20 K respectively. 



Comb ining IRAS and Spitzer observations of the Perseus molecular cloud, 



Schnee et al 



(120081 ) build a dust temperature map that shows the same range of values as within the 
TC. Such decomposition cannot be detected in our analysis as we may average it with the 
warmer envelope. We thus neither confirm nor question such a result at this time. Instead, 
we examine variations of dust temperature across the TC. 

The pixelwise dispersion in the observed emission ratios for the longer wavelengths 
may be accounted for by local variations of the dust properties and more particularly the 
BGs temperature. In order to analyze its spatial variations down to the smallest possible 
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angular scale, we build a map of the BGs temperature across the entire cloud using the 
temperature that best fits the ratio between the fluxes at 100 and 160 fim. The resulting 
map is shown on Figure O The mean value of the BGs temperature is 14.5 K and the 
spatial dispersion through the entire cloud is about ±1 K. With only two FIR data points 
at 100 and 160 fim we can not constrain the spatial variations of (3 as well and thus keep 
13 = 2 throughout this process. Had we assumed a lower value of /?, the absolute values of 
the temperature map would be larger but its structure would be the same. 

The BGs temperature dispersion appears to be related to the structure of the cloud. 
Cold filaments seem to be in the middle of a warmer medium. The temperature map is well 
correlated with the dust emission in the FIR. The cold structures are probed by the high 
fluxes at 160 /im, as shown on Figure [7l However, the FIR dust temperature we are using is 
an effective temperature. It corresponds to a line-of-sight average where the value at each 
depth into the cloud is weighted by the local emission. Since the FIR dust emission steeply 
increases with dust temperature, the effective temperature is biased towards the largest 
temperatures, those at the surface of the clouds. Smaller values of the dust temperature 



(5-. 12 K) have been derived from sub-mi 



( IStepnik et al. 



2003 



Schnee fc Goodman 



imiter observations of dark clouds within the TC 



2005 



Schnee et al. 



20071 ). This is probably due to 



the longer wavelength nature of the data. 



3.4. Dust opacity and extinction maps 



We use the dust temperature map to build FIR dust opacity maps. For every pixel 
of our dust temperature map we compute the Planck function and convolve it with the 
IRAS 100 /im or MIPS 160 /im transmission filter. We combine these values and the dust 
FIR emission at 100 or 160 /x r n to b uild maps of the dust opacity at 100 and 160 /xm 



(see Fig. 8(b)). 



Padoan et al. 



(I2OO2I ) have computed the structure function scaling of a 
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(a) (b) 
Fig. 8. — Comparison between (a) the visual extinction map from 



Padoan et al. 



(I2OO2I ) based 



on NIR observations and (6) the dust opacity at 160 fim deduced from our dust temperature 
map. 



2MAS S extinction map of the TC complex using the methods proposed by lCambresy et al 



( I2OO2I ). Not surprisingly, there is a strong correlation between their extinction map and 
our dust opacity maps (see Fig. [H]). From these correlations, we can compute the dust 
emissivity, assuming it follows a law in z/^. We find tioq/A^ = (3.65 ± 0.09) x 10^^ mag^^ 



and TiQo/A^ = (1.34 ± 0.04) x lO^^ mag-^ (see Fig. [9(^ 



Boulanger et al 



(119961 ) derived, from COBE and HI observations at high Galactic 
latitude, a diffuse dust emissivity that corresponds to 1.18 x 10~^ mag~^ and 0.46 x 
0~^ mag~^ at 100 an d 160 fim respectively using Ay/Nn = 0.53 x lO^^^cm^ 



( ISavage fc Mathis 



19791 ). There is a factor of about 3 between the cirrus measurements 



and those fro m the TC. Howe v er, if we assume the dust temperature is 17.5 K within 



the cirrus, as 



Boulanger et al. 



(119961 ) measured it on FIRAS observations, then their 



cirrus DIRBE colors from 100 to 240 /xm are better adjusted with an emissivity law of 
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z/^ ''3. If we use this law for the TC, we then find riooM^, = (2.58 ± 0.06) x 10 ^ mag ^ 
and T^p,n/A„ = (1.08 ± 0.04) x 10^^ mag~^, a factor 2.2-2.3 above the cirrus value 



from 



Boulanger et al. 



(Il996l ). Variations of the hydrogen column density to visual 
extinction conversion Ay/Nu could be the reason to explain the differences revealed by 
compari ng our measurements, deduced from direct comparison to visual extinction, to 



those of iBoulanger et al. 



(119961 ). who have been using HI observations. However, FUSE 



observations have sho wi i that this ratio d oes ri ot ch ange from d i ffuse to translucent clouds 



( iRachford et al 



2OO2I) 



del Burgo et al 



fl2003h and 



Kiss et al. 



mm used ISOPHOT 
observation to measure the dust emissivity at 200 /im, T20Q/Ay = 0.28 x 10"'^ mag~^ for 
lines of sight towards the TC, and T2oo/Ay = 0.2 — 1.0 x 10"'^ mag"'^ for other galactic lines 
of sight. The difference between these estimates and our values might only arise from the 
wavelength dependance of the FIR dust emissivity. It might also be due to the fact that 
their numbers are based on observations most often covering a too small area to provide a 
background outside of the cloud. An increase of the FIR emissivity by a factor of a few 
relative to the diffuse IS M has also been observed in the cold component of the Taurus 



Molecular Cloud TMC-2 fidel Burgo fc Laureijs 



2OO5I ). They interpreted this increase as an 



indication that the transiently heated VSGs are not present within the cold component but 
rather coagulate on bigger grains and that ice mantles form on dust. 



Figure 9(b) shows the variation of the visual extinction to FIR depth ratio at 160 /xm 



as a function of the BGs temperature for one IRIS plate that covers the TC. The BGs 
emissivity significantly increases with decreasing dust temperature. The correlation slightly 
varies over the entire complex, from one IRIS plate to another. This may arise from the 
difference in the dust temperature range that is spanned within each plate. The IRIS 



plate that is used for the correlation plot of Fig. 9(b) is the one that exhibits the largest 
contribution of pixel below 14 K. This correlation be tween low BGs temp erature and high 
FIR emissivity is in agreement with the results from 



Schnee et al. 



((20081) for the Perseus 
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Fig. 9. — Correlation plot of (a) the FIR depth at 160 fim as a function of visual extinction, 
(b) the visual extinction to FIR depth ratio at 160 fim as a function of the BGs temperature 
and (c) the dust emission at 160 fim as a function of visual extinction. 

Molecular Cloud. 

In order to confirm the good correlation between MIPS 160 /im data and the dense and 
cold structures in the cloud, we compare these FI R observations to the visual extinction 
map. On average, the resolution of the map from 



Padoan et al. 



(12002f ) is 3.4' but varies 



depending on the local number density of stars. We convolve the IR observations by a 



variable beam before the comparison is done. The correlation plot is shown on Figure 9(c 



The most extinguished structures are those with the most FIR emission and are the cold 
filaments on the dust temperature map. The good correlation between dust emission at 
160 fim, low BGs temperature and the extinction confirms that the dust FIR emission is a 



good tracer of dens e reg ions within mo 



Langer et al. 



fll989h and 



Abergel et al. 



ecula r clouds as it has previously been shown by 



199J). 



3.5. Dust modeling 



We compare the a verage r aid- t o far-infrared SED of the TC to that of the Diffuse 
Galactic Medium from iFlageyl (120071 ). At these wavelengths, the Taurus dust emission has 
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a steeper spectrum than in the diffuse galactic medium and this can be related to both 



different grain size distri 
using th e dust model of 
model of 



jution and tempe r ature . We explore both of these interpretations 



Desert et al. 



Compiegne et al 



(I2OO8I ). an updated version of the well-known 



(119901 ) to fit the average SED of the TC from 60 to 240 /im. The VSG 
and BG abundances and the interstellar radiation field (ISRF) intensity are free parameters 
in the model. Hereafter, we express the ISRF intensity in units of the Solar Neighborhood 
far-UV field: 1.6 x 10~^ erg.s"^.cm~^ integrated from 912 to 2000 A (Habing unit). The 
best fit model spectrum is shown as a solid line in Figure [51 

The corresponding ISRF is ~ 0.2, which is significantly low and may not represent 
the real mean incident radiation field. Indeed, the model does not take into account the 
observed increase in the far-IR dust emissivity. Therefore, the lower dust temperature, as 
compared to that of the diffuse Galactic medium, is entirely accounted for by reducing the 
radiation field intensity. However, for a given grain temperature the power radiated by dust 
in the far-IR is proportional to the dust emissivity while the absorbed power scales with the 
dust absorption coefficient in the UV/optical and the radiation field intensity. The balance 
between absorption and emission implies that for a larger dust emissivity, the same grain 
equilibrium temperature is obtained for a proportionally larger absorbed power. Then, 
assuming that the UV/optical dust absorption coefficient does not change, the previously 
stated ISRF intensity of ~ 0.2 would rather be in the range 0.5 to 0.6 if the model included 
the observed increase in far-IR dust emissivity. 

For each dust component (PAHs, VSGs and BGs), the model emission scales with 
the quantity of emitting grains (i.e. the dust mass in each component) and the grain 
heating power (i.e. the radiation field intensity). For a given emission to be fitted, a low 
value of the radiation field translates in a higher dust mass in each dust component. As 
a consequence, the dust total abundance is overestimated, though the relative fractions 
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of each dust component are we 
1.1 times that of 



Desert et al 



1 fitted. The abundance ratio of BGs and VSGs is about 



(119901) and 2.3 times that of 



FlagevI (120071 ) for the diffuse 



Galactic medium. However, this abundance ratio rehes mostly on the 60 /xm emission from 

lis difference in the 



Stepnik et al. 



(120031 ) 



VSGs. The increase of the FIR dust emissivity can be related to t 
VSGs to BGs relative abundance. Within a dark cloud of the TC, 
observed an increase of the sub-millimeter emissivity coupled with a decrease of the VSGs 
abundances, with factors relative to diffuse medium larger than ours, and interpreted them 
by introducing grain-grain coagulation into fluffy aggregrates as an important process inside 
the cloud. In the following section, we discuss evidence for spatial variations in PAHs and 
VSGs abundance. 



4. Local variations of MIR observations 

4.1. Selection of the sub region 

Zodiacal light is the dominant component of the diffuse emission at 8 and 24 /im and 
may contribute significantly to the observed spatial variations. We select a region, 2.6 
square degrees in extent, away from the brightest zodiacal bands, that exhibits filamentary 
structure at both 8 and 24 /xm (see Fig. [TOl) . At the shortest wavelengths (8 and 24 /^m), 
three filaments are clearly visible, forming a "three fingers hen-foot" shape. At longer 
wavelengths (60 /im and above), one of these filaments is significantly brighter than the 
two others, to the extent that it remains the one fully visible at 100 and 160 /xm. This "3 
fingers" structure is also visible at the same position within IRIS 12 and 25 /im images 
(see Fig. [T]and[2]). We thus assume it is part of the TC and not of the foreground zodiacal 
emission. 

We perform two cuts through the filaments. The positions of these cuts are shown in 
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Fig. 10. — Observations from 8 to 160 fim of the "3 fingers" sub-region at the IRAS 100 fim 
resolution of 4.3': (a) IRAC 8 fxm, (6) MIPS 24 /im, (c) IRIS 60 /im, (ti) MIPS 70 fxm, (e) 
IRIS 100 /im, and (/) MIPS 160 ^m. The dashed lines correspond to the two cuts made 
across the filamentary structures, which are shown on Figure [Til The units are in MJy/sr. 



Figure [TO] and their orientation is from high to low right ascension. Each map is at the 
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Fig. 11. — Relative MIR and FIR intensities along the two cuts as shown on Figure [TOl (a) 
is the cut #1 and (b) is the cut #2. Continuous lines are Spitzer values: (red) MIPS 160 
/xm, (green) MIPS 24 fim and (blue) IRAC 8 fim. Dashed lines are IRIS values: (red) IRIS 
100 fim and (green) IRIS 60 fim. Dotted lines are the temperature as deduced from the BG 
temperature map (Fig. [6]). The hatched areas represent the borders of each filament where 
the local minima have been measured to estimate the background. 

resolution of IRAS 100 /im (4.3') so we take a one pixel wide straight line for each cut. The 
variation of MIR and FIR emission along these cuts is shown on Figure [TTl On each cut, 
we identify several filaments thanks to the common presence of peak at every wavelength, 
from 8 to 160 fim. The relative position and amplitude of these peaks significantly change 
from one wavelength to another, from one peak to another. From east to west, we label 
filaments #1 and #2 on the first cut, filaments #3, #4, #5 and #6 on the second cut (see 



Fig. 10(a) to identify them on the IRAC 8 fim image). If we look back at the observations 
of this region, we see that filament #2 and filament 7^6 appears to be the same "finger" 
or spatial feature, while the "finger" or structure that contains filament #1 splits into two 
sub-features that correspond to both filament #3 and filament #4. Filament 7^5 seems to 
be part of a structure connected to filament #6, judging by the IRAC 8 fim observations. 
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Projection effects may significantly affect such associations. 

4.2. Average colors of the sub-region 

We measure the MIR to FIR colors of this field as we did previously on the entire 
cloud (see Section [3^ . At longer wavelengths (at 60 fim and above), the correlations are 
still good, except when MIPS 70 /xm is involved because of the significant noise of the data. 
At shorter wavelengths (IRAC 8 vs MIPS 24 and MIPS 24 vs IRIS 60), the dispersion is 
still significantly higher preventing a clean correlation and measurement of colors for those 
wavelengths. The measured FIR colors within this selected region are shown in Table [2] as 
well as those for the whole TC. The FIR colors of the "3 fingers" sub- region are globally in 
agreement with those from the entire TC. 

4.3. Small spatial scale variations in dust SEDs 

In order to constrain the variations of the dust properties on small spatial scales, we 
measure the colors of each filament independently. We first determine the borders of each 
of these structures on the cuts by looking for the minima values within the hatched areas 
(see Fig. [TTi) . For a given filament, the minima position may change from one wavelength 
to another by at most 2'. We then interpolate a straight line between the east and west 
borders to estimate the background that we remove from the filament intensity. We finally 
integrate the background subtracted intensity between the two borders of the filament. The 
final colors of the filaments, relative to MIPS 160 fim, are shown in Table H] and compared 
to the TC reference model in Fig. [121 For each filament, we give the background-subtracted 
brightness within the MIPS 160 channel as well as the corresponding visual extinction. 
The uncertainties on our measurements are dominated by the photometric uncertainties of 
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the observations: 10% on IRAC8 and MIPS24 and 15% on M 



( iReach et al 



2005 



Engelbracht et al. 



2007 



Gordon et al. 



2007 



PS160, IRIS60 an d IR: 



Stansberry et al. 



SlOO 



2003). We 



give these absolute photometric errors but they are not relevant for the relative analysis of 
small-scale variations. We rather draw the reader's attention on the relative errors, which 
are much smaller. For instance, the standard deviations at 8 and 24 fim, which include 
effects from structure and optical depth, are about 0.08 and 0.13 MJy/sr respectively, a 
factor 70 and 50 below the mean brightness at these wavelengths. 

The filaments colors clearly show that there is a difference between filaments #1, 
7^3 and #4 as a group and filaments #2 and #6 as another. The first three present 
IRIS100/MIPS160 and IRIS60/MIPS160 colors similar to those of the global TC (see Table 
E]). The last two present IR1S100/MIPS160 and IRIS60/MIPS160 a factor from 1.5 to 3 
above those average colors. At shorter wavelengths, the difference between these two sets of 
filaments increases as the IRAC8/MIPS160 and MIPS24/MIPS160 colors change by a factor 
from 2 up to 14. Filaments #2 and #6 present the highest colors at these wavelengths 
too. The changes we observe in the IR colors of these filaments occur within low opacity 
structures. The visual extinctions are at most a few magnitudes (see Table |2]). 

The difference at longer wavelengths can be related to different BGs temperatures, 
filaments #2 and #6 being "hotter" than the others, as shown in Figure [TT] There is 
a difference of about 0.5 K between the structure that corresponds to filaments #2 and 
#6 and the one that includes filaments #1 and #3. This is in agreement with what we 
expect from our previous analysis of the FIR brightness evolution as a function of dust 
temperature on the whole molecular cloud. Filament #5 appears to be in between the two 
groups of filaments but, as it is the faintest of our sample and as it presents a significant 
shift regarding the peak position, especially at 8 fim, we think it is not a conspicuous result 
and will not discuss it any further. 
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Fig. 12. — Spectral energy distribution of the filaments within the "3 fingers" sub- region of 
the TC, from 8 to 160 /xm, normalized at 100 iim. Reference model spectra for the Galactic 
diffuse medium and the TC are also plotted. 

5. Discussion 

The color variations observed across the TC may result from either a change in dust 
relative abundances or in the excitation efficiency. We explore these two possibihties by 
comparing the measured colors of the filaments to those of the dust model. Our aim is 
to understand whether the variations of IR colors within the filaments can be related to 
variations of only interstellar radiation field (ISRF) intensity x or if dust size distribution 
variations are required. In the following section, we use the average SED of the TC 
as a reference, both in terms of the ISRF intensity and dust size distribution, which is 
described by the relative abundances of PAHs, VSGs and BGs. As diagnostics, we use 
color-color diagrams that combine IRAC8/MIPS160, MIPS24/MIPS160, IRIS60/MIPS160 
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Compiegne et al. 



(120081 ) as a function 



Fig. 13. — Color-color diagrams from the model of 
of the ISRF intensity, VSGs and PAHs relative abundances, (a) is IRIS60/MIPS160 vs 
IRIS100/MIPS160, (b) is MIPS24/MIPS160 vs IRIS100/MIPS160, (c) is IRAC8/MIPS160 
vs IRIS100/MIPS160. The arrowed lines show the evolution of the colors as a function of 
increasing ISRF (from x = 0.3 to x = 2.8) for TC average dust abundances. The dotted 
lines show the same evolution for a factor 5 larger {red) or lower (green) relative abundance 
of the PAHs (squares) or VSGs (triangles) . Blue diamonds are the colors measured on the 
filaments. Blue stars on (a) are the average colors measured on the entire cloud and on the 
"three fingers" sub-area. The large blue crosses give the error-bars. 
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and IRIS100/MIPS160 ratios (see Fig. [13]). Colors predicted by the model as well as colors 
measured within the TC are plotted in these diagrams. 



5.1. Variations of the ISRF intensity 



We first attempt to interpret the observed color variations of the filaments by varying 
X from 0.3 and 2.8 times the reference ISRF. The dust size distribution is that of the 
average SED of the TC for the VSG and BG abundances (hereafter Yysofi and Ybg^)- 
The PAH abundance (hereafter Ipah.o) is arbitrarily fixed at YpAHfilYsGfl = 4.0 x 10^^ 
so that the IRAC8/MIPS160 ratio predicted by the model is within the range of observed 
IRAC8/MIPS160 r atios. As a comparison, the original PAH abundance in the model of 



Desert et al 



(Il990l ) is YpAH,o/yBG,o = 6.7 x 10^^. We show the model IR colors as a function 
of X with black arrowed line in the color-color diagrams. The IRIS100/MIPS160 color is 
the most sensitive to the variations of x- It varies by a factor ~ 2.5, while IRIS60/MIPS160 



varies by a factor ~ 1.4 (see Fig. 13(a)). MIPS24/MIPS160 and IRAC8/MIPS160 colors 



are almost independent of x (see Fig. 13(b) and 13(c)). 



The measured colors of the entire TC and the "3 fingers" sub-region are in agreement 
with those of the model. However, filament #1 is the only one whose colors match in every 
diagram, within uncertainties, the modeled colors for the TC average dust abundances. A 
few other filaments colors (e.g. IRIS60/IRIS100 and IRIS100/MIPS160 ratios of filaments 



#4 and 7^2 on Fig. 13(a) ) are in agreement with the model and only require different values 
of X- The X variations may be related to the extinction variations within the sub-region. 
However, the observed dispersion of the other filaments colors cannot be interpreted as due 
to variations of the ISRF intensity only. 



The variations of x may be responsible to modifications of the PAHs electric charge 
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that can be traced by IRAC8 color (iFlagey et al 



20061 ). Within the range of x used here, 



we have estimated the variations of the IRAC8/MIPS160 color due to changes of the PAHs 



charge to be 4% at most. Therefore, they are small relative to those observed in Fig. 13(c^ 
and we rule out this effect as a cause of the IRAC8/MIPS160 variations. 



5.2. Variations of the dust size distribution 

We change the dust size distribution by multiplying or dividing YpAHfi or YvsGfi by 
a factor 5. In order to disentangle between the effects of PAHs abundance variations and 
VSGs abundance variations, we modify only one abundance at a time {Ypah or Yysg) while 
the two others {YysG and Ybg or Ypah and Ybg ) are fixed. We thus create four different 
dust size distributions for which the abundance ratios relative to the reference model 
iYpAH/YpAHfi : Yvsg/Yvsga : Ybg/Ybg^) are (5:1:1), (1:5:1), (1/5:1:1) and (1:1/5:1). We 
compute the IR colors for each size distribution as a function of x and plot them in Figure 
[T3l All but one of the observed IR colors are within the range of values that the modeled 
colors span with the modified abundances. 

Not surprisingly, the IRIS100/MIPS160 color is almost independent of the dust size 
distribution variations. It is unaffected by a change of the PAHs abundance and is only 
slightly modified by a change of the VSGs abundance. The variations of the dust size 
distribution are more significantly traced by the other colors. The IRIS60/MIPS160 color 
appears to directly trace the variations of the VSGs abundance relative to BGs and is 



almost independent of the PAHs abundance variations (see Fig. 13(a)). However, the 
precise contribution of the VSGs to the IRIS 60 /im is still debatable. Besides, there may 
be different types of BGs with different equilibrium temperatures (see Section IX^ . Some 
of them could have a slightly higher equilibrium temperature than the single temperature 
that we assume in the model and thus could make a significant fraction of the IRIS 60 
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/iin emission. This statement does not apply to the emission in the MIPS 24 /im channel, 
which is only due to stochastically heated grains. In the TC, the contribution of the BGs 
to MIPS 24 fim emission is negligible. The MIPS24/MIPS160 color is mainly dependent on 
the VSGs abundance variations even though it is not completely independent of the PAHs 



abundance variations (see Fig. 13(b)). The last diagram shows that the IRAC8/MIPS160 



color is almost entirely dependent on the PAHs abundance variations (see Fig. 13(c)). We 
use these two last colors as tracers of the PAHs and VSGs abundance variations within the 
TC. An increase of the VSGs abundance is required for filaments #6, #2 and to a lesser 
extent #4; a decrease of the VSGs abundance is required for filament #3; no change of the 
VSGs abundance is required for filament #1 and 7^ 5. An increase of the PAHs abundance 
is required for filaments #2, 7^6, #5 and to a lesser extent #4; a decrease of the PAHs 
abundance is required for filament 7^3; no change of the PAHs abundance is required for 
filament #1. 

Thus, the colors of filaments #2, #3, #4 and #6 require an increase/decrease of both 
the PAHs and VSGs abundances while colors of filament #1 are in agreement with the 
reference size distribution. It is clear that one has to invoke PAHs and VSGs abundance 
variations by a factor of a few to explain the observed dispersion of the IR colors within the 
TC. In order to get quantitative measurements of the abundances variations, we compute 
the best fit of the parameters (x, Ypah, ^vsg) in Table [31 Even though the uncertainties on 
the radiation field intensity might be significantly higher than those on the abundances, we 
also give the best-fit values of x in Table [3l Relative to the TC reference, both the PAHs 
and VSGs abundances are lowered by a factor ~ 3-4 in filament #3 and increased by a 
factor ~ 2 — 5 in filament 7^2 and 7^6. The required increase is smaller within filament ^4 



(a factor ~ 1.2-2). In Figure 13(a), variations in the IRIS100/MIPS160 color corresponds 
to x/Xo ~ 1 for filaments #1 and ^3 and x/Xo ~ 3 for filaments #2 and 7^6. If this 
interpretation is right, we observe that the most illuminated filaments are those with the 
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largest proportion of PAHs and VSGs. Reciprocally, the less illuminated filaments are those 
with the smallest contribution of PAHs and VSGs. An alternative possibility - that would 
need to be tested with further modelling - is that variations in the IRIS100/MIPS160 color 
trace changes in the VSGs size distribution. We also note that the variations of the PAHs 
and VSGs abundances are not equal. Within the uncertainties, filaments #1, #3 and #6 
exhibit the same relative variations for both PAHs and VSGs abundances, while filaments 
7^2 and #4 require a relative increase of the PAHs abundance by 3 and 1.5 times that of 
the VSGs abundance, respectively. 

These abundances variations occur on sub-parsec distances as the characteristic width 
of the filaments within the "3 fingers" sub-region is about 20' or 0.8 pc (see Fig. [TTj) . They 
also occur within translucent sections of the cloud as the visual extinction of the filaments 
is about a few magnitudes at most (see Table Sj). Since the filaments were selected in an 
unbiased way, except for avoiding the zodiacal light, we could infer that similar dust size 
distribution variations take place also through the entire TC. 



5.3. Evolutionary processes 



The dust size distribution variations observed within the TC can be placed into a broader 
context. For instance, 



Boulanger et al. 



([19901) and 



Bernard et al 



19931) observed enhanced 



Miville-Deschenes et al 



f l2002h 



abundances of the smallest dust particles at cloud surfaces, 
detected an increase of the PAH abundance within a specific H I velocity component in 
Ursa Major. However, these previous studies, based on IRAS and ISOCAM observations, 
lacked the sensitivity and field of view of the recent Spitzer observations. The IRAC 8 /im 
and MIPS 24 fim data provide much higher spatial information on the PAHs and VSGs 
emission than in past studies. On the one hand, they reveal new structures within the cloud 
that remain unnoticed either in the Ay or MIPS 160 /xm map. On the other hand, within 
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the "clean" ecliptic latitude strip that we use for this analysis, there is no IRAC 8 fim and 
MIPS 24 ij,m counterpart to the bright structures in the Ay and MIPS 160 /xm images. 

We relate the observed abundances variations to mass exchange between the small 

dust particles and the biggest grains. In the average model of the TC, the PAHs and VSGs 

contribute to 9% of the total dust mass. Within filaments ^^3 and #6, which represent the 

two extreme cases of our analysis, the contribution of PAHs and VSGs to the total mass of 

dust ranges from 3 to 27%. The mass fraction of PAHs and VSGs is 12% in the original 
I 1 

model of iDesert et al.l (119901). We compute the equivalent fraction of carbon mass that is 



involved in those exchanges. We appl y a carbon-silicate dust mass proportion of 37:63 which 



1984 



Weingartner fc Draine 



200l|). We 



fits the diffuse ISM extinction curve (jPraine fc Lee 
use a carbon abundance [C/H]gas = 161 ± 17 and [C/H]totai ~ 330 C atoms per million 
H determii ied along translu cent lines of sight, even though these values are still subject 



to debate (ISofia et al 



2004 



y, and references therein). For the TC, the average fraction of 
carbon within the smallest dust particles (PAHs and VSGs) is about 26% relative to the 
total mass of carbon in dust particles and 13% relative to the total mass of carbon in both 
dust and gas. Within filament ^^3, where the relative abundance of PAHs and VSGs is 
the lowest, these figures are about 8% and 4% respectively. Within filament 7^6, where 
the contribution of PAHs and VSGs is the largest, the values are about 72% and 37% 
respectively. Therefore, a significant fraction of the dust and carbon mass moves in and 
out of the narrow size range of stochastically heated PAHs and VSGs. These dust/carbon 
mass exchanges occur at the cloud surface, for moderate values of the visual extinction, well 
below the threshold where ice formation occurs (e.g. extinction thresholds of 4.3 ± 1.0 mag 



for CO^, 3.2 ±0.1 rnag fo r H2O and 6.7 ± 1.6 mag for CO in the Taurus complex. 



Whittet 



2003 



Whittet et al. 



20071). 



MathisI (|l99J) and 



Whittet et al 



(I2OO4J ) suggested that changes in the UV extinction 
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curve could be accounted for by the formation of hydrocarbon coatings on grains small 

enough (size < 20 nm) to contribute to the 2200 A bump and far-UV rise of the extinction 

curve. Such a process could also account for the low abundance of PAHs and VSGs within 

the 160 fim bright regions of the TC. The timescale for coagulation depends not only on 

the density but also on the relative velocities between g rains. Grain relat ive velocities are 

I 1 

set by turbulent motions and depend on the grain size. lYan et all (|2004| ) have calculated 
the relative grain motions arising from magnetohydrodynamic turbulence as a function 
of ISM physical conditions. They find that the largest velocities are always those of 
BGs. In that case, small grains are expected to condense on large grains faster than they 



mutually coagulate (jPraine fc Anderson 



19851 ). However, small grains can also be efficiently 



decoupled from gas motions by large velocity gr a dients at sm all scales, known to exist 



in turbulence and not included in the 



Pety &: Falgarone 



YanetaL 



(120041 ) work (IFalgarone fc Puget 



1995 



20031 ). If small grains acquire relative velocities comparable to those of 
BGs, PAHs and VSGs coagulation would occur faster on small grains than on BGs alone, 
because the total surface in small grains is larger. 

Shattering of hydrocarbon grain mantles in grain-grain colli sions is a likely process to 



expla in the local enhancements of PAHs and VSGs abundances ( Serra Diaz-Cano fc Jones 



19961). Here, 



20081 ) . Shattering occurs above some collision energy threshold ( iJones et al. 
turbulence may play a signicant role in creating the necessary relative velocities between 
grains. The correlation between dust emission at 8 fim and 24 fim and specific structures in 
the gas velocity thus appears like a promising perspective for a future study. Once specific 
measurements of velocity gradients are available, the spatial scales over which the observed 
variations occur could be converted into their equivalent timescales. 
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6. Conclusions 

We have presented Spitzer images of the Taurus Complex. We have combined these 
images with IRAS and DIRBE observations to characterize the diffuse emission across the 
cloud from mid to far-IR wavelengths. The sensitivity, spatial resolution and wavelength 
coverage of the observations permits analysis of small spatial scale variations of the dust 
emission. The data interpretation highlights evidence of dust evolution occurring within 
the translucent sections of an archetype example of quiescent molecular clouds. 

The comparison between IRAS and MIPS observations at 60, 70, 100 and 160 //m 
reveals a strong spatial correlation between wavelengths although we can identify some 
variations. The average colors of the TC provide us with a FIR spectral energy distribution 
that is steeper than that of the Galactic diffuse medium. This difference suggests colder 
BGs. The average dust temperature is about 14. 5K with a dispersion of ±1K across the 
cloud. The FIR dust opacity is tightly correlated with the extinction derived from 2MASS 
stellar colors. The FIR opacity to extinction ratio we infer from this correlation is a factor 
2.2 larger than the average value for the diffuse ISM. This increase in opacity and the 
attenuation of the radiation field both contribute to the lower emission temperature of the 
large grains. 

The structure of the TC diffuse emission significantly changes at MIR wavelengths. 
Due to zodiacal light contamination, the Spitzer maps only offer a partial view at the MIR 
emission from stochastically heated PAHs and VSGs. Outside the zodiacal bands, where 
the contamination is small, the 8 and 24 /^m maps are similar but there are no 8 and 24 /xm 
counterparts to the brightest 160 /xm emission features. The 8 and 24 /xm images reveal 
filamentary structure that is strikingly inconspicuous in the 160 /im and extinction maps. 
The analysis of these filamentary structures shows significant variations within the MIR to 
FIR colors, from one filament to another and within one filament, over sub-parsec distances. 
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We present model calculations that quantify the dependence of MIR to FIR colors on 
the radiation field intensity and the abundance of stochastically heated particles. To match 
the range in observed colors, we have to invoke variations by a factor of a few of both the 
interstellar radiation field intensity and the abundance of PAHs and VSGs. 

Our results should motivate theoretical studies of the evolutionary processes, grain 
shattering and coagulation, which could account for the observed changes in the spectral 
energy distribution of the dust emission. The correlation between dust emission within 
IRAC 8 fim and MIPS 24 /xm and specific structures in the gas velocity would allow 
determination of the potential role of turbulence in driving changes in the abundance of the 
smallest dust particles. 

The authors would like to thank Laurent Cambresy for providing the 2MASS extinction 
map. They also thank Paul Goldsmith, Vincent Guillet and the anonymous referee for their 
questions and comments, which have significantly improved the paper. 

This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. 
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Table 1: Average DIRBE colors and DIRBE 100 /xm to IRIS 100 /xm ratio across the TC 

Ratio Value 

DIRBE60/DIRBE100 0.17 ±0.02 

DIRBE100/DIRBE140 0.29 ± 0.02 

DIRBE140 /D1RBE240 0.81 ± 0.05 

DIRBElOO/IRISlOO 0.82 ± 0.01 



Table 2: MIR-to-FIR colors of the "three fingers" region and the whole TC. 



Color 


Value 






"Three fingers" 


TC 


IRIS60/MIPS70 


0.63 ±0.22 


0.46 ± 0.07 


IRIS60/IRIS100 


0.18 ±0.02 


0.20 ± 0.02 


MIPSTO/IRISm 


0.31 ±0.15 


0.35 ±0.05 


I Ri Sim/ MI p Sim 


0.27 ±0.01 


0.25 ±0.01 
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Table 3: Dust components abundances and interstellar radiation field intensity used within 
our model for the best fit of the filaments. 



Filament 


#1 


#2 


#3 


#4 


#5 


#6 


Yvsa/yvsGfi 


0.92 ± 0.05 


1.8 ±0.2 


0.25 ± 0.04 


1.3 ±0.1 


0.48 ±0.25 


3.3 ±0.1 


Ypah/Ypah,o 


0.92 ±0.14 


5.2 ±0.8 


0.33 ± 0.04 


1.9 ±0.3 


3.3 ±0.5 


3.8 ±0.6 


X/Xo 


0.9 


2.5 


0.8 


1.1 


1.1 


3.2 



Tabic 4: MlR-to-FlR colors of the filaments within the "three fingers" region, normalized 
to MIPS 160 /xm, in xlO~^. The background subtracted peak flux within MIPS 160 iim 
channel as well as the equivalent visual extinction are also given. 



Fflament #1 #2 #3 #4 #5 #6 

IRAC8/MIPS160 9.1 49 3.4 18 34 35 

M1PS24/M1PS160 8.7 32 2.6 14 8.5 32 

IRIS60/MIPS160 51 77 30 58 82 120 

IRIS100/MIPS160 290 490 270 320 310 500 

MIPS160 (MJy/sr) 37 6 66 12 4 9 

Av (mag) 2.1 0.3 3.8 0.7 0.2 0.5 



